NEUTRON MEASUREMENT METHOD FOR DETERMINING POROSITY OF 
A FORMATION SURROUNDING A BOREHOLE 



BACKGROUND OF THE INVENTION 

5 This invention relates to the determination of a characteristic of earth 

formations penetrated by a borehole, and more particularly to the determination of 
formation porosity through neutron measurements. 

In hydrocarbon exploration and production, it is of prime importance to 
determine (a) if a given earth formation contains hydrocarbon, (b) the amount of 

10 hydrocarbon within the formation, and (c) the producibility of the hydrocarbon in 

place within the formation. The amount of hydrocarbon present within a formation 
is a function of the pore space or the "porosity" of the formation. In drilling wells 
for the production of hydrocarbons, it is desirable to measure the porosity of each 
prospective hydrocarbon producing formation penetrated by the borehole. It is even 

15 more desirable, for economic and prospective reasons to determine the porosity of 

prospective formations during the actual drilling of the borehole. 

Over the past decades, many technologies have been used to measure or 
estimate formation porosity from a borehole. One of these technologies is based on 
a system that contains an isotopic source that emits fast neutrons, and an axially 

20 spaced detector that responds to the flux of impinging thermal neutrons resulting 

from the interaction of fast neutrons with nuclei within the borehole and formation 
in the vicinity of the borehole. The basic concept of this system is predicated on the 
fact that (a) hydrogen is the most effective moderator of fast neutrons because of its 
low atomic weight, and (b) most hydrogen found in earth formations is contained in 

25 liquid in the pore space of the formation, either as water or as liquid hydrocarbon or 

gas. The detector is axially spaced from the neutron source such that for a given 
borehole condition, the count rate recorded by the thermal neutron detector 
decreases as the volumetric concentration of hydrogen, or porosity increases. 

Dual detector neutron porosity systems have been introduced to minimize 
30 the effects of the borehole upon the measurement of formation porosity. US 
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3,483,376 and US 5,767,510 disclose two thermal neutron detectors that are spaced 
axially at different distances from the source of fast neutrons. The ratio of the 
responses of the two detectors varies with formation porosity, yet is less sensitive to 
borehole parameters than the count rate from either two individual detectors. The 
5 ratio is therefore the measured parameter used to compute porosity. Historically, 

this ratio has been formed from the response of the detector closest to the source, or 
the "near'' detector, divided by the response of the detector farthest from the source, 
or the "far" detector. 

However, accuracy problems remain with these two detectors systems since 
10 the response of said neutron porosity tool significantly varies with the density of the 

formation to be logged. Furthermore, this response is also a function of the thermal 
capture cross section (sigma). Minimizing the sigma response of a neutron porosity 
measurement is usually achieved by enclosing the neutron detectors in a highly 
absorbing thermal neutron shield through which only epithermal neutrons may 
1 5 penetrate. In that way, only epithermal neutrons are detected, resulting in very little 

sigma response. 

Careful positioning of detectors with respect to the source can minimize 
density effects. It has been shown that a minimum in density response occurs at a 
certain unique distance from the source depending on the source energy (see Scott, 
20 H.D., et .ah, 1994, "Response of a Multidetector Pulsed Neutron Porosity Tool", 

paper J, in 35^^ Annual Logging Symposium Transactions of the Society of 
Professional Well Log Analysts). Detectors placed at the point of minimum density 
sensitivity have been shown to have very little density response. 

However, this technique has various drawbacks, among which is a strong 
25 restriction on the placement of the detectors, which in turns, leads to important 

mechanical constraints. Furthermore, the far detector still displays significant 
density sensitivity, 

SUMMARY OF THE INVENTION 

It is thus an object of the invention to propose a method for determining 
30 porosity of a formation surrounding a borehole providing design flexibility while 

greatly improving the dynamic range of the neutron porosity measurement; 
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minimizing environmental corrections and eliminating differences in shale response 
among various tool designs. 

In order to achieve this goal, a neutron measurement method is proposed for 
determining porosity of an earth formation surrounding a borehole comprising: 

5 - conveying a tool along said borehole, wherein said tool comprises a 

source of neutron radiation and at least one detector axially spaced from 
said source; 

generating measured detector response for said at least one detector that is 
indicative of neutron radiation from said source interacting with said earth 
1 0 formations; 

operating said measured detector response with a predetermined 
mathematical equation and thereby obtaining corrected detector response 
that is independent of the density of said earth formation; and 

determining porosity of the earth formation surrounding the borehole from 
1 5 said corrected detector response. 

Preferably, the predetermined mathematical equation comprises multiplying 
the measured detector response by a correction factor that depends on the density of 
the formation. Advantageously, the measured and corrected near detector responses 
comprise a near detector count rate, said measured and corrected far detector 
20 responses comprise a far detector count rate. 

Advantageously, said mathematical equation is of the form: 

wherein CRcon- is the corrected detector response, CR is the measured detector 
response, P is the detector sensitivity to density and p is the formation density. 

25 It is also an object of the invention to propose a system for determining 

porosity of an earth formation surrounding a borehole comprising: 

(a) a borehole tool comprising a source of neutron radiation and at least one 
detector; and 
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(b) a computer for computing measured response of said detector thereby 
obtaining a measure of the porosity of the earth formation surrounding the 
borehole, whereby: 

- said measured response from said at least one detector is indicative of nuclear 
5 radiation from said source interacting with said earth formation; 

- said measured response of said detector is operated with a predetermined 
mathematical equation using said computer to obtain corrected detector 
response that is independent of the density of the formation; said corrected 
detector response being indicative of the porosity of the earth formation 

1 0 surrounding the borehole. 

The advantage of these method and system is a more accurate neutron 
porosity measurement having desirable response characteristics. Moreover, this 
neutron measurement method applies to neutron porosity measurements made with 
various means of conveyance: wireline (WL), logging-while-drilling (LWD), or 
1 5 logging-while-tripping (LWT). 

BRIEF DESCRIPTION OF THE DRAWINGS 

Additional objects and advantages of the invention will become apparent to 
those skilled in the art upon reference to the detailed description taken in 

20 conjunction with the provided figures: 

figure 1 shows a neutron porosity measurement tool; 

figures 2a and 2b represent uncorrected and corrected count rates 
according to the method of the invention; 

figures 3a and 3b represent uncorrected and corrected near to far detector 
25 ratios according to the method of the invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Figure 1 illustrates a detector neutron porosity device embodied as a 
logging- while-drilling tool. In this example of a system according to the invention, a 
source of fast neutrons 1 , a near detector 2 and a far detector 3 are positioned within 
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a drill collar 4, which will be referred to as the LWD tool. In another example of a 
method and system according to the invention, the tool can utilize more than two 
axially spaced detectors, thereby allowing single or pairs of detectors responses to 
be selected to optimize the limits of the prospected zone. The LWD tool 4 is 
suspended by means of a drill string 5 within a borehole 6 penetrating an earth 
formation 7. The upper end of the drill string 5 is suspended at the surface of the 
earth 70 preferably from conventional rotary drilling draw works (not shown). A 
drill bit 8 is affixed to the lower end of the LWD tool. The drill string 5 is rotated by 
means of a kelly 9 thereby rotating the LWD tool 4 and drill bit 8, and thereby 
extending the borehole downwardly as a result of the cutting action of the drill bit. 
A preferably conventional drilling fluid system is employed to remove cuttings 
formed by the rotating drill bit, and to maintain hydrostatic pressure within the 
borehole. 

In the example of figure 1, the near detector 2 is axially spaced at distance d 
from the neutron source 1 and the far detector 3 is axially spaced at distance D from 
said neutron source. The neutron source, near and far detectors are pressure sealed, 
preferably within the wall of the tool 4, thereby isolating these elements from the 
borehole environment, and also allowing for a preferably coaxial channel within the 
tool 4 through which the drilling fluid flows. The drilling fluid exits the drill bit 8 in 
a manner well knovm in the art. The neutron source is preferably an isotopic source 
that emits fast neutrons. This could be a chemical source such as a mixture of 
Americium and Beryllium (Am-Be, 4.5 MeV), or alternately, Californium-252 
(^^^Cf). The neutron source could also be a neutron generator using a d-d or d-T 
reaction. The pulsed neutron generator consists of a high voltage supply (typically 
70 to 100 kV) and a vacuum tube in which deuterium and tritium ions are first 
generated and then accelerated onto a target containing tritium and deuterium. The 
reaction between a deuteron and a tritium nucleus results in the emission of a 
neutron with an energy of about 14 MeV. 

The near detector and the far detector are preferably sensitive only to very 
low energy neutrons, that is "thermal" or "epithermal" neutrons, Helium-3 detectors 
are responsive to both thermal and epithermal neutrons, while a Helium-3 detector 
wrapped with a layer of cadmium is responsive primarily to epithermal neutrons as 
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is well known in the art. In the example where a ratio between the far detector and 
the near detector responses is implemented, it is preferred that the far detector 3 be 
more sensitive to thermal neutrons for statistical reasons, since the flux of thermal 
neutron at the position of the far detector will be considerably less than the thermal 
neutron flux at the near detector. 

Due to the method of the invention, when a far and a near detectors are 
provided, the relative position of the far and near detectors can be varied with 
respect to the neutron source, which is very convenient since it allows mechanical 
adaptability. Referring to figure 1, the near and far spaced detectors can both be 
positioned above the neutron source at preferably axial spacings, respectively. 
Alternately, either the near or far spaced detector can be positioned above the 
neutron source, and the other detector (or both) can be positioned below the neutron 
source. 

Power supplies and control and data conditioning circuitry for the detectors 
are contained preferably within the LWD tool 4. The counting rates of the near 
detector 2 and far detector 3 are preferably telemetred to the surface of the earth 70 
by means of a mud pulse telemetry system, not shown, or other suitable telemetry 
system known in the LWD and MWD (Measurement While Drilling) art. 
Altemately, detector responses or the count rate data can be recorded and stored 
within a memory means (not shown), preferably located within the LWD tool 4, for 
subsequent retrieval when the LWD tool is returned to the surface of the earth 70, 
using a computer (not shown) to display the near and far count rates, subsequently 
creating a log of the porosity of the formation as a function of depth within the 
borehole. 

In general, detector responses depend on both hydrogen index & density of 
the formation surrounding the borehole (assuming no sigma effects, where sigma is 
the neutron macroscopic capture cross section of the formation, see e.g. The book 
"Well logging for Earth Scientists", by Darwin Ellis (Elsevier)). 

Then the count rate (CR) for a given detector can be expressed by: 

CR^e-^"^ xe-PP eqn(l) 
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Where: HI = formation hydrogen index, i.e. the hydrogen concentration of 

the formation relative to that of water. 

p = formation density 
a = detector sensitivity to HI 
5 P = detector sensitivity to density 

and an exponential dependence on HI and density is assumed. The sensitivities a 
and P are constants for a given detector spacing and neutron source energy (in 
particular, P varies depending on chemical or neutron generator sources but with the 
appropriate coefficient, this differences can be overcome to attain accurate porosity 

10 measurements). The sensitivity P is determined for a given detector by observing 

the count rate response while varying the formation density but keeping the 
formation HI constant. The sensitivity a is obtained by doing the reverse operation. 

Then the equation : CRcorr = CR x e^^ e"""^ 

removes the density effect and isolates the HI response. This equation can then be 
1 5 applied to at least a measured detector response or both to a near and a far detector 

count rates. 

Thus, according to the invention, corrected detector count rates (Ncrcoit and 
FcRcorr rcspcctively) in the example wherein two detectors (near and far) are 
implemented can be obtained by: 

20 NcRcorr = NcR X C^"^ ~ C"""^""' 

FcRcorr = FcR X t^'^^'^ ~ e-f^^"' cqii (2) 

wherein Ncr and Fqr are the measured near and far detector count rates 
respectively. 

According to the method of the invention it is possible to use a single 
25 measured detector count rate and correcting it with said method in order to obtain, 

with an appropriate processing method known in the art, the formation porosity. It is 
also possible to use a pair of near and far detectors and then to use either both 

corrected near and far detector count rates or only one corrected count rate and a 
second measured count rate in order to obtain a count rate ratio. However, using 
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both corrected near and far detector count rates will be preferred in order to obtain 
the most accurate formation porosity. 

In the example where both near and far detector count rates are corrected, 
this will lead to the corrected ratio: R = — ^Aj^^hl^ This ratio is relatively insensitive 

to density and will respond only to hydrogen index. In particular, as represented on 
figures 2a and 2b, readings in shales should not yield exaggerated porosities caused 
by density effects. That is, the apparent porosity of a 45 pu shale using the 
uncorrected ratio (fig.2a) is about 82 pu, whereas the density corrected ratio (fig.2b) 
yields the correct porosity of 45 pu. 

To further explain, figure 2a shows a far detector measured response prior to 
density correction. The deviations of the sandstone and dolomite points from the 
limestone porosity response are due primarily to density effects (especially for 
dolomite). The alumina point, representing a dense shale response and well 
separated from the limestone response line, displays even more density effects. In 
Figure 2b the far detector count rates have been corrected for density effects, 
according to the calculation method of the invention. As a result, the dolomite and 
alumina points lie on the limestone response line, demonstrating that density effects 
have been removed .The sandstone points move closer to the limestone response 
line but do not overlay the limestone response due to residual effects having to do 
with fast neutron transport. 

As represented in figures 3a and 3b, the effects of density correction are 
even more apparent in terms of the ratio of near to far count rates. The uncorrected 
Ncr/Fcr ratio (fig.3a) shows significant density effects for sandstone, dolomite and 
the alumina point compared with the limestone response (that is, the sandstone, 
dolomite and shale points lie well away from the limestone porosity response line), 

N . 

On the other hand, the density-corrected ratio R = '''''''' is insensitive to density 

effects and will respond only to hydrogen index HI. As a result, the dolomite and 
shales points lie on the limestone response line (the sandstone responses are much 
closer to the limestone response line but exhibit residual effects). 
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Applying the density correction to the detector count rates has several 
benefits. First, the dolomite and alumina points lie on the limestone response line 
and thus will need no lithology correction, that is, the porosity response in 
limestone, dolomites and shales is very nearly identical. Secondly, the sandstone 
points lie closer to the limestone response line and hence will require a smaller 
lithology correction to place them on the reference limestone response line. Thirdly, 
and most importantly, the dynamic range of the measurements has been 
dramatically improved, especially at the higher porosities, by removing the density 
effects. That is to say, the dynamic range (0-^100 pu) is about 5 for the uncorrected 
ratio and about 20 for the density-corrected ratio. Better dynamic range results in a 
more accurate and statistically precise porosity measurement. 

The method of the invention can be applied to an existing LWD tool design 
to "tune" the porosity response of the tool so that it exhibits a specified amount of 
density effect. The amount of density effect exhibited by a given design is 
dependent on the source/detector spacing and can be quantified by the beta 
coefficient of eqn (1). By adjusting the value of beta in the density correction of eqn 
(2), density sensitivity can be increased, decreased or removed altogether to attain a 
specific density response. The ability to tune the porosity response of the tool in this 
way makes it possible to deliver a response that is independent of the tool design. 
Thus, a consistent and uniform porosity response can be obtained from tools having 
intrinsically different responses. 

The method of the invention has been described in reference to an LWD tool 
but provides equally effective density correction in wireline systems. 
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